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We introduce a model for directed percolation with a long-range temporal diffusion, while the 
spatial diffusion is kept short ranged. In an interpretation of directed percolation as an epidemic 
process, this non-Markovian modification can be understood as incubation times, which are dis- 
tributed accordingly to a Levy distribution. We argue that the best approach to find the effective 
action for this problem is through a generalization of the Cardy-Sugar method, adding the non- 
Markovian features into the geometrical properties of the lattice. We formulate a field theory for 
this problem and renormalize it up to one loop in a perturbative expansion. We solve the various 
technical difficulties that the integrations possess by means of an asymptotic analysis of the diver- 
gences. We show the absence of field renormalization at one-loop order, and we argue that this 
would be the case to all orders in perturbation theory. Consequently, in addition to the characteris- 
tic scaling relations of directed percolation, we find a scaling relation valid for the critical exponents 
of this theory. In this universality class, the critical exponents vary continuously with the Levy 
parameter. 
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I. INTRODUCTION 

Directed percolation (DP) H Q is the generic 
model for nonequilibrium systems which exhibit a con- 
tinuous phase transition into a unique absorbing state. 
DP describes the temporal-directed spreading of a non- 
conserved agent in a certain medium. The agent might 
be bacteria in a case of epidemics in populations, fire in 
a burning forest, or water in a porous rock. The spread- 
ing phenomenon is characterized by two competing pro- 
cesses: relying on the medium conditions, the agent may 
multiply itself or decay at a constant rate. Depending 
on the balance between these two processes, the spread- 
ing may continue forever or die out after certain time. If 
the agent is not allowed to appear spontaneously, in the 
latter case the system is trapped in the absorbing state, 
a configuration where the stochastic fluctuations cease 
entirely and the system cannot escape from. These two 
regimes of survival and extinction are typically separated 
by a continuous phase transition characterized by the DP 
critical exponents. 

The DP universality class is extremely robust, as a 
whole range of theoretical models seems to belong to it. 
Some examples include heterogeneous catalysis Q , chem- 
ical reactions , interface depinning '3| , the onset 
of spatiotemporal chaos 0, flowing sand (Ifl] and self- 
organized criticality The robustness of the model 
has led to the conjecture that two-state spreading pro- 
cesses with short-range interactions generically belong to 
the DP class, provided that quenched randomness, un- 
conventional symmetries, and large scales due to memory 
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effects are absent 0, . 

In the context of epidemics, DP describes infection pro- 
cesses without immunization and where the disease is 
only transmitted to nearest neighbors by direct contact. 
This can be understood realizing the problem of DP on a 
{d + l)-dimensional lattice, where each lattice site is con- 
sidered as an individual which can be infected (active) 
with probability p or healthy (inactive) with probability 
1 — p. An infected individual recovers at the next time 
step with probability 1 and is ready to be reinfected with 
the same constant probability p. The susceptibility to 
infection is then independent of previous infections, and 
this ensures the absence of immunization. 

Therefore, in order to make a realistic description of 
epidemics the effect of immunization as well as long-range 
interactions should be taken into account as modifica- 
tions of the DP model. Immunization can be added to 
DP by considering a probability for subsequent infections 
different from the first infection probability This 
non-Markovian feature changes the universality class 
of the model to the one corresponding to dynamical 
percolation, also known as a general epidemic process 
(GEP) J^, 16] . The phase diagram of this model displays 
a curve phase transition line connecting the GEP and DP 
critical points Along this line the same universal- 

ity class as GEP is found. A horizontal phase transition 
line also separates the GEP phase from a supercritical 
DP behavior. The absence of scaling along this line has 
been shown for the case of 1 -I- 1 dimensions and later 
generalized to + 1 dimensions [T^ l20l | . 

On the other hand, an epidemic model with long-range 
interactions was first suggested by MoUison 21]. This 
model was studied as a generalization of DP which in- 
cludes spatial long-range interactions where the spread- 
ing distances follow a power-law distribution given by. 
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II. THE MODEL 



V{r)dr ^ r 



'dr. 



(1) 



where d is the spatial dimension and cr is a control pa- 
rameter. Asymptotically, as r ^ oo, Eq. (Q) is equivalent 
to a Levy distribution and a is called the Levy expo- 
nent. In this sense we can say that the particles perform 
Levy flights j^^- The claim that the critical exponents 
should vary continuously with a |23j| was confirmed by 
theoretical renormalization-group analysis [24l| as well as 
by extensive numerical simulations |25l l2q . More re- 
cently these results have been generalized to branching- 
annihilating Levy flights and to the pair annihilation 
reaction A + A ^ [2^ . So far all these studies have 
assumed dynamic processes which are local in time. In 
the case of epidemics, it is assumed that the infection 
happens instantaneously in time. 

In order to make one step forward in making a more 
realistic model for epidemics, in this paper we set out 
to study a non-Markovian modification of the DP prob- 
lem which includes waiting times, or incubation times, 
between the infection and actual outbreak of the disease 
in a population. We assume that these incubation times 
T are distributed asymptotically as r — > oo 



!F{T)dT oc 



1 



dr. 



(2) 



Here the Levy parameter k > is a free parameter that 
controls the characteristic shape of the distribution. For 
r — > we assume J^{t) is a smooth function of r. For 
simplicity we also assume that the dynamic processes are 
local in space, which means that the infection can only 
spread by contact with nearest neighbors. 

The characterization of the universality class of this 
problem has remained an open problem in the literature 
since it was first suggested in a previous work 26] in 
1999. This is mainly because of the technical difficul- 
ties that arise in the field-theoretical description when 
long-range waiting times are introduced. In the present 
work, we derive a field theory for this problem and cal- 
culate the critical exponents by means of systematic per- 
turbation theory and the e expansion. This paper is or- 
ganized as follows. In Sec. ^ we propose a convenient 
approach to derive the field-theoretical action, through 
a generalization of the method introduced by Cardy and 
Sugar in 29]. We dedicate Sec. IIIII to study and an- 
alyze the mean-field predictions of our theory. Subse- 
quently in Sec. IIVI fluctuation effects are taken into ac- 
count via renormalization-group methods. The various 
difficulties that emerge through the rcnormalization pro- 
cess are managed by studying the asymptotic behaviour 
of the integrals involved in the renormalization process. 
Finally, in Sec. we write the renormalization group 
equations and compute the critical exponents at one loop. 
In Sec. IVII devoted to the conclusion, we argue that our 
results are valid to any loop order. 



A. Master Equation 

In order to derive the field-theoretical action for the 
problem of DP with incubation times, we first consider 
the master equation formalism. Directed percolation can 
be interpreted as a reaction-diffusion process of identical 
particles in a d-dimensional lattice, where multiple oc- 
cupation is allowed. We call P{a,t) the probability that 
the system will be at a given microstate a. The dynamics 
of such a system usually is described by a master equa- 
tion governing the temporal evolution of the probability 
distribution P{a,t), which its general form is given by 



dP{a,t) 
di 



The system goes from the microstate a to the microstate 
(3 with a constant transition rate Ra^fj. A naive gen- 
eralization of this to processes involving transitions with 
incubation times would be 



'^^^ dt'Rp^^{t - t')P{P,t') - 

V / dt"R^^p{t" ~t)P{aA), (4) 
^ Jt<t" 



where the transition rates Ra-^p are time-dependent 
functions. But Eq. Q is wrong, as the probabilities 
P(/3, t') do not refer to mutually exclusive events for dif- 
ferent times t' . Indeed, Eq. Q does not conserve the to- 
tal probability P[oi, t) . In fact Eq. describes the 
dynamics of particles, which disappear from the lattice 
at time t', until they reappear at time t > t' . This does 
not correspond to the dynamics with incubation times 
that we are trying to model. 

To write a correct master equation one should add to 
the right-hand side of Eq. |0J) an infinite number of terms 
which will take into account the nondisjoint nature of the 
events, so the master equation is replaced by an infinite 
hierarchy of coupled equations for multitime joint prob- 
abilities. 

In order to avoid dealing with the difficulties inherent 
to a master equation formalism, we propose to adopt an 
alternative way to find the field-theoretical action. We 
will generalize a method first introduced by Cardy and 
Sugar in order to show that directed bond percolation 
was in the same universality class of Reggeon field the- 
ory m. 
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FIG. 1; Example on a (l+l)-climensional lattice of possible 
temporal long-range infection vectors from a lattice site {x, t) 
to nearest-neighbor sites in space. 

B. Action of DP with incubation times 

In order to provide a model for epidemics with long 
incubation times, we consider our system on a (d -I- 1)- 
dimensional lattice. We represent the spreading of the 
infection on the lattice through vectors (see Fig. ^| . An 
infection vector between a lattice site (x, t) and another 
site {x',t') is present with probability p{x' — — t), 
with t < t'. The temporal coordinate t indicates the 
preferred direction, and therefore the orientation of the 
infection vectors is always in the direction of the increas- 
ing time. The vectors can only connect nearest neigh- 
bors in space, but their range in time depends on the 
incubation times distributed as Eq. Q . Considering this 
model, the problem of epidemics with long incubation 
times can be interpreted as a temporal long-range di- 
rected percolation problem. We propose now to write 
a field theory for this model, through a generalization 
of the Cardy-Sugar method. We define the connectivity 
function G{x2,t2] xi,ti) as the probability for two given 
lattice sites {xi,ti) and {x2,t2), with ti < t2, to be con- 
nected irrespective of the other sites. We say that two 
sites are connected if there is an infection vector present 
between them. Following Cardy and Sugar G can 
be written as 

G{x2^t2\xi,ti) — TT:a{x2,t2) ]^ [l+p{;x—x,t' 

links,t' >t 

— t)a{x' ,t')a{x,t)]a{xi,ti). (5) 

The commuting operators a{x,t) and a{x,t) act on each 
site {x, t) of the lattice, and their algebra is defined as 



(J? — ia, (J? = za, (6) 



Tra(a;,i) = Tra(a;,<) = 0, (7) 



Tr[a(a;,f)a(x,t)] = 1. (8) 

We remark that Eq. ||SJ) is identical to the one obtained 
by Cardy and Sugar in 2!||, except for the fact that in 
our case the probability p is not a constant. 

The physical features of the problem and the details of 
the dynamics are included in the effective lattice deter- 
mined by the infection vectors. We define a matrix V ^ 
which will contain this information as follows: 

[1 + — t' — f)a(x', t')a(x, f)] 

linhs >t 

= exp(^^a(a;',t')F(x' -a;,i' -t)o(a;,t) j. (9) 

^ x,t x' .t' ' 

Therefore the matrix Y will contain the information of 
the temporal long-range processes. In order to complete 
the generalization of the Cardy-Sugar method for the 
problem of DP with incubation times, we assume that 
Y can be decomposed into a short-range part Ysix,^) 
and other part Yiix, t) which will be long range in time, 

V{x,t) = Vsix,t) + Viix,t). (10) 

Yi (x, t) contains the long-range temporal dependence and 
a factor with a spatial dependence, which is short ranged. 
Therefore, we can assume that the leading behavior of 
Yi{x,t) is as follows: 

y^M-^^ (11) 

with a proportionality factor that is some short-range 
function of x. We consider an expansion of the Fourier- 
Laplace transform of Yg in a small momentum k and a 
small energy E, 

Vs{k,E) = E(i - ^(^^)' -Et+- ■ • V.(x,t), (12) 

x,t \ 

where 

Vs{k,E)=c-ciE~C2k^ + -- - . (13) 

In the case of long-range temporal processes considered 
here, the moment (t) is divergent and we cannot perform 
an expansion in E. Instead, we compute the Laplace 
transform of JF(t) as 

/■OO 

/ e^'^^T{t)dt E" + const + regular terms. (14) 
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Therefore, the Fourier-Laplace transform of the long- 
range contribution Vi{x,t) will involve a E'^ dependence 
multiplied by the Fourier expansion of a spatial short- 
range factor, 

Vi{k,E) - +const + • •■)(! - 6fc2 -!-•••)• (15) 

We should notice that this is valid for values of k > 0. 
Keeping the most relevant terms in a small-fc and -E 
expansion, the Fourier-Laplace transform of V is given 
by 

V{k, E) = c[l - riE - rafc^ - rE'' + Oik^E'')]. (16) 

Applying Gaussian integrations in Eq. @, we can be 
written as 

exp ^ ^ a{x', t')V{x' -x,t' - t)a{x, t) ) 

^ x,t x',t' ' 

/oo 
d(/)#exp[^(a;',i')(^)^V(a;,t) -a0-a0], 

(17) 



III. MEAN-FIELD APPROXIMATION 

A. Critical exponent /3 

In this section we will find the mean-field values of the 
critical exponents of the theory. If fluctuations effects 
are neglected, the field can be interpreted as a density 
field, and consequently above criticality it scales as 

0(x|pc-p|'^, V>Vc- (21) 

We start by finding the classical equations of motion. In 
order to do so we consider a variation of the action in 
Eq. H20(l . with respect to the fields and 0. If we define 
the Lagrangian density as 

H = ^(af + T9t - i?oV2 + ro)0+ ^u.M'' - ^uo'^''/', (22) 
the corresponding equations of motion are 

d'i<i> + M + ^wo('/'' - 2#) - i^oV^f^ = 0, (23) 



where the operator V ^ \s. given by 

V'^ = cr\\ -\- ridt - rzV^ + rd^ 



(18) 



Finally using Eqs. © and (|17|l . and replacing the result 
in Eq. we obtain 



G = TTa{x2,t2)a{xi,ti) 

D(j>D^exp(- / dtd'^x[^{XV)-^(l) - a4> - acj)] 



(19) 



Performing the trace operation in Eq. (|19(l and after ap- 
plying a rescaling of the fields, we find the effective action 
of the problem of directed percolation with incubation 
times. 



S = / dtd'^x 



(20) 



Here rp ccpc— p and the fields and cf) depend on space 
and time. For k = 1, Eq. H20|l is reduced to the DP 
case. Thus, in our problem we consider < k < 1. The 
generalization of the Cardy-Sugar formalism, in contrast 
with the master equation method, does not simply add 
the rates for any possible infection from time t to another 
time t' in the future, but it properly takes into account 
the multiple counting of events which may occur if there 
is more than one way of reaching a given state at time t' . 



roi + -uo{24,cj) - 02) + {-l)-d^4> - DoV^ = 0, (24) 

after using Eq. H22|) . A solution where 4> — would be 
equivalent to not considering the noise fluctuations in 
the Langevin-like equation. From Eq. H24|) we see that 
^ = is indeed a classical solution, as long as Eq. is 
verifled: 



(25) 



A particular solution of this equation of motion can be 
obtained if we neglect the temporal and spatial depen- 
dence of the field 4> — that is, a mean-field approximation. 
Therefore, Eq. (|25|) becomes 



1 2 



(26) 



giving two solutions, for tq > (below criticality), 

= 0, p<pc, (27) 
and for tq < (above criticality), 
, -2ro 



Mo 



P> Pc 



(28) 



Therefore, in a mean-field approximation (j) tq, and 
from here we obtain that 



P^'^ = 1. 



(29) 
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Critical exponents i'± and 



In order to calculate the exponents v± and i^y, we will 
analyze the scaling behavior of the correlation function 
G^^'^\x,t), around the Gaussian fixed point when the 
interaction terms in the action are neglected. Below crit- 
icality, we do not expect any longer a temporal exponen- 
tial decay of G*^^'^-* (a;, t), as happens in the case of pure 
DP, but a power-law behavior. This is due to the fact 
that the infections can happen at at very large times. 
However, the spatial decay of G*^^'"'^^ (x, t) is an exponen- 
tial decay, since the spatial diffusion is short ranged. 
Therefore G^^'^-* (a;, i) decays exponentially in the limit 
of large x and as a power law in the limit of large times. 
We proceed to write the Fourier-Laplace transform of 
G(i^i)(a;,i) as 



7+200 



dE 



r,Et 



2ni S« + Dok^ 



If we make the change of variables E — E'r^^'^ and k 

k 



(30) 



1 / 1/2 

-ttsTq I the correlation fmiction can be rewritten as 
follows: 



1/2 



(31) 



-1/k 



and 



Consequently, time and space scale as t 

— 1/2 

X ^ , respectively. At criticality ro — 0, and there- 
fore any temporal and spatial scale is divergent. We can 
then define the critical exponents ly^, which describes 
how the spatial correlation length diverges at criticality, 
the exponent , describing the divergent behavior of the 
temporal correlation length, and the dynamic exponent 
z = vnjv^^ such that 



MF 



,MF 



1 

1 

K 



2 

K 



(32) 



The value of these exponents are given at a mean-field 
level, as we have derived them neglecting the interactions 
in the action in order to compute G*^^'^-* (x, t). 

Next, we will find how G'^'^-' (x, t) decays below criti- 
cality. We should notice that the Laplace transform in- 
volved in Eq. H30|l cannot be solved exactly. Thus, in Ap- 
pendix^ we compute how this integral behaves asymp- 
totically in the limit of t — > oo, finding 



G(i'i)(x,t) 



1 



^1+^ 



{P < Pc)- 



(33) 



We should compare this result with the DP problem, 
where G^^'^-* {x, t) decays exponentially below pc- At crit- 
icality we can perform similar calculations setting a; — > oo 
(see Appendix 0, and obtain 



G(i'i'(oo,t) 



(P =Pc)- 



(34) 



Consequently we find that G^^'^\x^t) behaves with dif- 
ferent power laws at criticality and below criticality. At 
criticality, G*-^'^' follows a power-law decay given by an 
exponent: 



6^'^ = l-K. 



(35) 



In this way we have shown that at a mean-field level, 
the critical behavior can be described by continuously 
varying exponents with the Levy parameter k. For k = 1, 
we recover the DP exponents. 



IV. FIELD-THEORETICAL ANALYSIS 



In this section we will include the effect of fluctua- 
tions in our analysis, and therefore a mean- field approach 
cannot be considered any longer. Instead we will apply 
field-theoretical techniques which will allow us to perform 
the calculation of the critical exponents below the upper 
critical dimension. We start by computing the canonical 
dimensions for the various quantities appearing in the 
action in Eq. H20|l . simply by considering the dimension- 
less nature of the action. In addition, the time-reversal 
symmetry (^ — > — ^; (j) — > —(j)), still valid in this problem, 
suggests the use of equal canonical dimensions for both 
fields, 6 and 6. Therefore, 



(36) 



and the dimensions of the fields depend on the Levy pa- 
rameter K. The canonical dimensions of the diffusion 
constant and the coupling constant uq are 



and 



"0 



= tj(3K-l)/2^-d/2^ 



(37) 



(38) 



respectively. In order to calculate the upper critical di- 
mension dc^ we express the canonical dimension of the 
coupling constant in terms of momentum units only, as 
follows 



ul 



D, 



(3k-1)/k 




^ ^(6K-2)/K-d 



(39) 



Hence, we see from Eq. that the coupling constant 
becomes dimensionless at the value of the upper critical 
dimension dc, 



dr. = 



6k - 2 



(40) 



below which the fluctuation effects become important. 
We should notice that Eq. H40|l gives a negative dc when 
K < 1 /3. For these values of k, a mean-field theory rather 
than a field-theoretical approach should be implemented. 
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Vertices 





FIG. 2: Propagator and vertices for the problem of directed percolation with incubation times. 



Consequently, in this section we only consider 1/3 < k < 
1. 

The Feynman rules for the propagator and the vertices 
of the theory are formulated in Fig. El The propagator 
G^'^'^^k, E) is represented by a straight line, and its ex- 
pression can be obtained from the free action, taking the 
Laplace transform of time, and the Fourier transform of 
the spatial dimensions into momentum space, given as 



1 



£;« + Dok^ + ro 



(41) 



We have neglected in Eq. (|41() the linear term tE with 
respect to E'^, since in the low-energy limit {E 0), 
the latter term is dominant. The main difference with 
respect to DP is the modification of the propagator due 
to the long-range temporal infections, given by the non- 
Mar kovian term E'^. Notice that the vertex interactions 
are not altered with respect to DP. 

In what remains of this section we proceed with the 
renormalization of the theory. We will apply mass, field, 
and diffusion constant renormalizations to absorb the di- 
vergences of the two-point vertex function F^^'-'^'. The 
divergences of F^^'^^ will be considered in the coupling 
constant renormalization, and finally we will renormalize 
the composite operator 



A. Mass and field renormalizations 

We start by writing the two-point vertex function F*^^'^^ 
at one loop. Figure 01 shows the diagrammatic expansion 
up to one loop — that is. 



dE' 



{2iTi) J {2itY 



[E'^ + Dok'^ + ro][iE - E'Y + Da{q - kf + ro] ' 

(42) 

This vertex function has two kinds of divergences. One 
kind may happen at d = 4 — 2/k, and we assume it 



is absorbed into a redefinition of the bare mass tq to a 
renormalized mass rn. A second kind of divergence may 
happen at d — 6 — 2/k, and it will be absorbed into a 
renormalization constant of the fields. 

We will work in a Laplace-Fourier space constituted 
by an energy E, considered positive and real, and a mo- 
mentum k. We define the normalization point (NP) such 
that the external energy is evaluated at an arbitrary scale 
E = while the external momentum is set to zero — that 
is, q — 0. We define the first renormalization condition 



(1,1) 



d{E'') 



= 1. 



(43) 



NP 



The renormalization of the fields defines the renormal- 
ization constants and Zr, such that 



(44) 



Nevertheless, due to the time-reversal symmetry, we can 
choose = Z^. The two-point vertex function, cal- 
culated by cutting off the external propagators to the 
correlation function G'-^'^-', is then 



-(1,1) 



r.(i,i) 



(45) 



Inserting this into the renormalization condition, 
Eq. (|43|l . we obtain the expression for the field renor- 
malization constant, 



<l> 



9F(i 



9(^;«) 



(46) 



NP 



Unfortunately, we were unable to evaluate the integral in- 
volved in the expression of F(i'i) in Eq. ^ exactly, and 
consequently we must rely on an analysis of its asymp- 
totic behavior at the singular points. Applying standard 
complex variable theory, we can see that the integral 
in Eq. H42|l presents two logarithmic branch points: if 
we write E'^ — e"''^^ , we identify one branch point in 
E' = 0. In the same way it is possible to see that E' — E 
is the second branch point. There are no poles in the first 
Riemann sheet, and therefore we consider the branch-cut 
topology shown in Fig. ^ 
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E, q 



r(i'i) = < + 

E"" + Doq^ + tq 





E\ k 




Wo / 

\ 






E - 


-E',q- 
E', k 


-k 








E 


-E\ q 


- k 



< 



FIG. 3: Feynman diagrams which contribute to the expansions of the propagator G*-^'^-* and the two-point vertex function 
r'^'^-* up to one loop. 



branch cut 



branch cut 



- integration contour 

FIG. 4: Branch points and branch-cut topology for the inte- 
gral in Eq. lUSll . 



In order to solve the integral /, we notice that there are 
two branch points present, one at u = and another at 
u = 1. We were unable to evaluate this integral exactly, 
and consequently in Apoendix lB ll we analyze the asymp- 
totic behavior of the integral at the possible points where 
divergences may occur. We find that there is no other 
divergence present in Eq. H49|l . except for the one reab- 
sorbed in the definition of the renormalized mass rji. The 
direct consequence of this is that the field renormalization 
coefficient remains constant — that is, = 1 + const — 
and under a suitable rescaling of the fields it is possible 
to redefine such that at one loop order we have 



Z,, = 1. 



(50) 



In conclusion, according to one-loop calculations, the 
field renormalization is not required in the theory and 



We will replace the bare mass by the renormalized 
mass Tji, assuming that the mass renormalization is al- 
ready done. For simplicity, when evaluating the integrals 
in Eq. (|46|l . we will set r/j = and 9 = 0. We start 
by performing the momentum integration which can be 
solved exactly in a straightforward manner, giving 



r(i 



with 



NP 



CSC — 
2 2 \ 2 



1 



d/2 



£,(l-e/2)Kj 



r/2-i 



{2TTi) 



(1 — u)'^ — u'^ 



E=c 
(47) 



(48) 



In Eq. g7l, e = 4 - and 5*^ = 2Tr'^/^ /T{d/2). We have 
made the change of variable E' — Eu, where E is real 
and positive. Taking the derivative with respect to E'^ 
and evaluating in the NP we obtain the expression for 
the field renormalization constant Zw,: 



ZV' = 1 



ul f dir 



D. 



d/2 



1 - 



C^"'/^/. (49) 



-(1,1) 



NP 



(51) 



NP 



which proves that the bare propagator is the full propa- 
gator for our theory. 



B. Diffusion constant renormalization 

In the case of the diffusion constant renormalization we 
proceed in a similar manner to that described in the pre- 
vious subsection, since in this case we also have the tech- 
nical difficulty that the integrals involved in the renor- 
malization cannot be calculated exactly. Then, we ana- 
lyze the asymptotic behavior of the integrals at the dif- 
ferent singularities in order to determine the divergences. 

We impose the renormalization condition 



(1,1) 



— Dr — ZdDq, 



(52) 



NP 
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FIG. 5: Perturbative expansions for the correlation function G*-^'^' and the corresponding vertex function F*-^'^' up to one 
loop. 



which defines the renormalized diffusion constant Dr and 
the renormalization constant Z^: 



Zt 



(53) 



NP 



We simplify our calculations setting from the beginning 
in our analysis rji — 0. After performing the momentum 
integration we obtain 



9r( 



1,1) 



Do 



Sd 



2 i{2TTyD^/^' 



-7rcsc( — ) 



1 



'B2 



(54) 



where the integrals Iei and Ie2 are given in Appendix lB 21 
by Eqs. (|B6|I and HB7|I . respectively. In Appendix IB 21 
we obtain the divergences of these integrals through an 
analysis of the their asymptotic behavior at the singular 
points. At the upper critical dimension, we find loga- 
rithmic divergences. The integrals cannot be solved an- 
alytically, and therefore we calculate the coefficients of 
the divergences. The results are given by Eqs. (jB14|) 
and IjBlSp . Inserting these results into Eq. H54|) . the ex- 
pression for Dn is given by 



Sd 



j^d/2-l 2K(d/2-3)+5(-27r)d 



xcsci —]^ [Fi{k) + F2{k)1 

2 / K€ 



where the functions Fi{k) and F2{k) are 

1 \ sin(7rK) 



Fi(k) = 2 2- 



K J sin"^(f k) 



0, (55) 



(56) 



and 



K sin (-1^) 



4k2 - 3k -f 1 
3k - 1 



37r 



(57) 



Although the approximations we have used to compute 
the behavior of the integrals are rather drastic, they lead 
one to obtain the coefficients of the divergences in a trans- 
parent way. 



C. Coupling constant renormalization 



Turning now to the coupling constant renormalization, 
we wish to construct and follow an equivalent procedure 
to tackle the integrations that will appear in this part 
of the renormalization process. We start by defining the 
renormalized coupling constant ur as 



ur 



-(24) 



-Z 



3/2p(2,l) 



NP 



-(2,1) 



NP 



(58) 



NP 



We should notice that in Eq. (|58|l we have used the fact 
that there is no field renormalization and therefore sub- 
stituted ^0 = 1 . In Fig. [S] we show the Feynman dia- 
grams contributing to r^^^^^) up to one loop. 

The dressed vertex function F'^'^^ truncated at two- 
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loop order is, then, 



dE' r d'^k 



{2m) J {2TTf 



1 



{E'- + Dak^)[{E^ - E'Y + Doiqi - 
1 

^ [{El +E2- E'Y + Da{qi + 92 - kf] 

+ perniutation(i?i, qi ^ E2, q2)- (59) 

The second Feynman diagram contribution, represented 
by a permutation(_Bi, qi <-> £'2, (72), is an integration with 
value equal to the first integration in Eq. H59|l . but with 
the external energy Ei and momentum qi exchanged 
with E2, and (72, respectively. In order to make the com- 
putation of the integral in Eq. H59|) easier, we choose the 
renormalization point NR as Qi = 52 — 0, E2 — 0, and 
El = Q. This choice is arbitrary and should not mod- 
ify the final results. Thus, evaluating Eq. (|59|) in the 
normalization point, we have 



^9 n ^ /■ dE' f d'^i 

°i 27rz 7 27r 



d'^k 

r 



{E"^ + i:>ofc2)[(£; - E'Y + Dok^Y 



(60) 



The factor of 2 comes from counting the contribution of 
the second Feynman diagram. We proceed by integrating 
first over the momentum, and we find 



p(2,l) ^ 



Uq n csc{dTr/ 2) Sd 



D, 



d/2 



{2ttY 



Iei 



(61) 



Remarkable, the Iei is the same integral which first ap- 
peared in the Dq renormalization. We then make use of 
the result in Eq. HB14|) and inserting it into Eq. (|59|l we 
finally obtain 



ul csc{d7T/2)Sd „ , X-'^'/' 



2K(d/2-3)+2(-27r)'i 



(62) 

We are now able to define and calculate the dimension- 
less renormalized coupling constant gn. Accordingly to 
Eq. ((23), it can be defined as 



9R 



D 



d/i' 



(63) 



where /i is a momentum scale and therefore related to 
C hy fi — CY^"^ jD^J^ . As we have used an energy scale 
in the renormalizations, it is convenient to write in 
terms of C, as follows: 



9R 



Ur 



D 



(d-e)/4 
R 



eK/4 



(64) 



Inserting the expansions of Dr [Eq. (|55|l ] and ur 
[Eq. (inU], into Eq. JHH) we have finally 



9R 



where 



Uq 



{d-e)/4' 



Uq 



d/4: 



eK/4 



Sd csc(d7r/2) 
2K(d/2-3)+7(27r)'^' 



-F,{^) 



F3{k) = (32 - dc)Fi{K) + dcF2{n) 



0, (65) 



(66) 



The next subsection is dedicated to study the final renor- 
malization step: the composite operator {(pcj)) renormal- 
ization. 



D. Composite operator renormalization 



We start by defining the renormalized composite oper- 
ator {(t>(t))R as follows: 



(67) 



Then the renormalized two-point correlation function 
with the insertion of the composite operator can be writ- 
ten as 



G^,^'^^^' = {{Yd^)RMR) 



(68) 



and the corresponding renormalized vertex function 
r^'^'^-* is defined by cutting off the external propagator 
in Eq. 

^ R 



We impose the normalization condition 



.(i,i;i) 
R 



1, 



(69) 



(70) 



NP 



and using Eq. I|69|l we obtain the expression for the renor- 
malized constant Z-jj,, 



(71) 



NP 



The next step is computing the unrenormalized vertex 
fimction F'-^'^'^-' evaluated on the renormalization point 
NP, chosen as q = and E = C^. The Feynman diagrams 
corresponding to the dress vertex function F^^'^'^-* up to 
one loop are shown in Fig. 

The simplest way to calculate F'-^'^'"'^' is, first, shift- 
ing above criticality where the renormalized mass tr is 
different from its value at criticality rRc- We define a 
parameter Aq = tr — trc as a measure of the departure 
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— <- 



r(i'i;i) 



< X < + 



q -^^^ 


E', k 












11/ — nj ^ 




Jz/ — JZ/ 


q — k 




q — k 




E', k 






/ — ^ 




<S( 




Vf.,.. 






^E-E' 


q — k 




q — k 



FIG. 6: Expansions of the correlation function g'^'^'^' and r'^'^'^' at one loop. The insertion of the composite operator 
is represented by a cross. 



from criticality, and we take the derivative of r'^-'^'^^ with 
respect to Ao, as follows: 



NP 



dAo 



NP 

dE' 



1 



X 



(27ri) {E - E'Y ~ E'' 
d<^k (-1) 

1 



{2TT)d [{E - E'Y + Dok^ + Ao] 



(72) 



Second, in order to evaluate r'^^^'^^ at criticality we take 
the limit Aq — > 0. After performing the momentum in- 
tegrations, we have 



where 



Sd 



NP 



8(2^) 



^(d-2Vcsc ^ ]I2, (73) 



d7r\ 



B/2+JOO £.,K(<i/2-2) £./)K(<i/2-2) 



E/2-ioo 



{2lTi) 



{E - E'Y - E" 



(74) 

In Appendix lB 3l we obtain the asymptotic behavior of I2 
at the logarithmic divergence [Eq. ljB19|l ]. Substituting 
this result into the expression oiT^^.U^) in Eq. we 
have finally 



Z7^ = l- 



sm 



Sd 



2K(d/2-3)+3(-2^)d 



{d-2) 



CSC — 

K€ V 2 



0. (75) 



In this way we complete the renormalization procedures 
required to absorb any possible divergent term in the 
vertex functions up to one loop. In the following section 
we will write down the renormalization-group equations 
and calculate the critical exponents. 



V. CALLAN-SYMANZIK EQUATION 



Having performed all the renormalizations required 
for the theory, we are now able to calculate the 
renormalization-group equation for F^'^' and r^'^'""^' at 
criticality and derive the critical exponents. We will de- 
rive the Callan-Symanzik equations considering a nor- 
malization scale fj, in units of momentum and make the 
change to ( through the relation fi-^ = §C^- Subse- 
quently the fi dependence will disappear when we express 
it in terms of physical quantities, such as energy and mo- 
mentum. We start by using the fact that the unrenor- 
malized F*^^'^' does not depend on the normalization scale 
introduced by the normalization point NP, having then 



dfi 



no, Do 



0, (76) 



and replacing now the total derivative with partial deriva- 
tives, we find 



R 



dgp. 



(77) 



where the flow equations, the beta and gamma functions, 
are define as follows: 



P{9r) = M 



uo,Do 



70 = M 



J no, Do 



ID^ fi 



dlnZo 



Juo,Do 

(78) 

The beta function can be calculated using Eq Ht)5|) . and 
it reads as 

R(n \ -'^ r( ^3r\ _ e ( Up ^„,^/4 
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FIG. 7: 13 function Pign) for k = 0.7 and e = 0.01. 



9r 



In order to write Eq. (|85|l in terms of the physical quanti- 
ties of momentum k and energy i?, we replace the deriva- 
tives on /i and Dji using the following identities: 



^'Tr dk 



and 



(89) 



which are easy to derive from direct application of the 
rule of chain. In this way, we eliminate the /i dependence 
in the Callan-Symanzik equation, rewriting it as 



where 



—Sd csc{dn/2) 

2K(<i/2-3)+7(-27r)<i 



(80) 



is a positive and finite function of k in the domain of 
interest — that is, ^ < k < 1. Up to first order in uq, 
can be written as 



Uo 



-£k/4 







and inserting it into Eq H65() we obtain 



^0 



(81) 



(82) 



Using this result we rewrite Eq. H79() obtaining an expres- 
sion of the beta function in terms of gj^: 



PidR) - 

which vanishes at the fixed point 



^B. + -9r + Oig%), 



9r 



(83) 



(84) 



The fixed point g*^ is an infrared-stable fixed point as is 
possible to see from Fig. The renormalization-group 
equation H77|) . evaluated at ^^ij, can be written as 



T'^^'^\^l,DR,UR)=Q, (85) 



where = l<j>{9*R) = and 7^ = 1d{9*r) at one loop is 
2 



Id 



[F^{K)+F^{K)]e + 0{t 



A solution of Eq. (|85|l is given by 



.(14) _ 



= DRfl'<fl-,— 



k 



E 



(86) 



(87) 



4-<- 



, + ,i) + i(2 + ,i)BA 



.(1.1) 



{k,E)^Q, 
(90) 



and applying standard methods to solve it; we obtain 

E 



T^^''\KE) 



fc(2+7D)/« 



(91) 



From this equation we can derive that E ~ fc^^+To)/", 
and considering the definition of the dynamic exponent 
z, E ^ t^^ ~ fc^, we obtain that 



Td 



which at one loop gives 

_ 2 ^ 2[Fi{k) + F2{k)] 

K kF^^k) 



(92) 



(93) 



This expression is valid at one-loop order in the pertur- 
bative expansion. The next step now is to determine the 
anomalous dimension. In order to do so, we calculate the 
two-point correlation function for large times < — > 00: 



G(i'i)(x,oo) 



d'^k 



jfc(2+7L)/K 
fc2-7;+7L 



(94) 



It is straightforward to find the value of the anomalous 
dimension through a simple comparison of the two last 
lines in Eq. finding that ??± = f + - 1). 

Using the result obtained for z in Eq. 1)92(1 and making 
7^ = 0, we have finally 



??± = z{l - k) + 2. 



(95) 



This new relationship at one-loop order between 77^ and 
z is a direct consequence of the absence of field renormal- 
ization. We want to compute now how the density of ac- 
tive particles behaves as x — > 00 — that is, the correlation 
function G^'^'(oo,t). Therefore, we turn our attention 
to Eq. where we can see that k ^ E^/^'^+^*o\ and 

as /c ^ 0, the vertex function Fp'^^ scales as 



^[K/(2+7|,)](2-7;+7i,)^ 



(96) 
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It is then easy to obtain the temporal dependence of Making use of the identities 



G^^'^^ at /fc = 0, using the behavior of the vertex function 



given by Eq. 

which gives a power law decay with the exponent: 



6=1- 



2 + 1*0' 



(97) 
(98) 



At one loop, because of the absence of field renormaliza- 
tion in the theory, 7^ = and therefore the value of this 
exponent coincides at one loop with its mean-field value 
found in Eq. H35|) . This is an expected result, in the sense 
that when field renormalization is not required, the bare 
propagator, valid to describe the density of active par- 
ticles at a mean-field level, is itself the full propagator 
of the theory when fluctuations effects are taken into ac- 
count. 



A. P, v^_, and i^u exponents 



In this subsection we shall investigate the 
renormalization-group equation for the vertex function 
p(i,i;i) at criticality. The results of these calculations 
will let us derive the critical exponents /3, and u^. 
The starting point again is the independence of the 
unrenormalized r^^'^'^^ on the momentum scale /Lt, which 
lets us write 



1t.(1,1:1) 



0. 



(99) 



We should notice that the parameter which accounts for 
the shift of criticality Aq = r/j — r^c is taken equal zero, 
and we will only use it at some point in the calculations in 
order to do dimensional analysis. The Callan-Symanzik 
equation for F^^'^'^^ reads as follows: 

(4+7^,-7;+7L^«^)r^^'^^^^=0, (100) 

where we have already evaluated the gamma function 



dlnZr 



have 



7i^ 



at the fixed point. At one loop 



kF^^k) tan(-|K) 



(101) 



We will maintain explicitly 7^ in the equations and make 
it zero in the end, with the only purpose of pointing out 
the direct consequences of the absence of field renormal- 
ization at one loop. Through dimensional analysis, we 
infer a solution of Eq. H100|l as follows 



.(1,1;1) 
R 



/k 
$ _ 



E 



(102) 



d 



d 



1 d 

\^dE 



and 



dk 

n ^ - V ^ 



(103) 



(104) 



we can replace the derivative on /i and Dji in terms of 
derivatives in momentum k and energy E, to obtain the 
Callan-Symanzik equation at criticality: 

A solution of this equation is given by 



E 



(106) 



Now, we can use the scaling form of the vertex 
function F^'^'* above criticality — that is, F^'"'^'' ~ 
k^-ti+^o f{k~^/''^ Ao, E-^/^ii Ao)— to obtain in an alter- 
native way the scaling form of rl^'^'^-*- 



.(1,1;1) 



R 

(1,1) 
R 



Ao=0 



(107) 



This equation in comparison with Eq. (|106|l allows us to 

find the value of the exponent iy\ as a function of 77 , 

4>4> 

and 7|), 

^ (108) 



2-7i, + 7l, 



At one-loop order this equation gives 
1 /8(2k-1) 



\Ktan(^K) 



(109) 

In addition, through the definition of z, we know that 
= zv^, and therefore we find 



kI 2 



Id 



■ 7L + ih 



which at one loop gives 



1 



+ 



8(2k- 1) 



K k'^ F3{k) ta,ii{^ k) 



(110) 



(111) 



We calculate now the /? exponent above criticality, by 
writing down how G^^'^-' {x, t) behaves in the limit of large 
times: 



df^k duj 



k 



Ao^-'A, 



27r 

lAoP''. (112) 
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In this way, above criticality we obtain the same relation- 
ship vahd for DP — that is, 



(d + ry^ -2). 



(113) 



We can write d — dc — e, and thus the value of the /3 
exponent at one loop can be calculated using Eq. (|113|l . 
given 



: 8(2k- 1) , , , , kF^(k) 



tan(7rK/2) 



2 4 

,2\ 



+ 0{e'). (114) 



Nevertheless, in our theory we have in addition an extra 
relationship given by Eq. (|95|) because of the absence of 
field renormalization in the theory. Using both Eqs. 1)95(1 
and H113() . we obtain, at one- loop order, 



2/3 = v\\{l - K) + dv^. 



(115) 



The existence of this relationship makes the exponent 
/? dependent on the value of and and therefore 
reduces the independent critical exponents from three to 
two independent critical exponents, with respect to the 
DP theory. Therefore, we find a new scaling relation at 
one loop for the problem of DP with incubation times. 
In the following section we argue that this result is valid 
to any loop order in perturbation theory. 



VI. DISCUSSION 



In this paper we have formulated and solved a field 
theory for a non-Markovian model of directed percola- 
tion with the inclusion of long-range temporal diffusion, 
which in a context of epidemics can be interpreted as 
incubation times. The incubation times are distributed 
following a Levy distribution, while the spatial diffusion 
as well as the interactions are short ranged. We first draw 
the attention to the fact that the conventional approach 
of writing a master equation, in order to apply later a 
second-quantized formalism, is not a convenient way to 
find the field-theoretical action. This is mainly due to the 
fact that the master equation has an infinite number of 
terms, a consequence of the nondisjoint nature of events 
for different times. Instead, we have proposed an exten- 
sion of a method introduced by Cardy and Sugar [2^. 
where we included the details of the long-range temporal 
diffusion in the effective lattice determined by the infec- 
tion vectors. Following this approach, in a rather simple 
way we have found the action of the problem. 

Second, we found already at a mean-field level that 
the critical exponents vary continuously with the Levy 
parameter, signaling the existence of a new universality 
class. We also found at a mean-field level that the two- 
point correlation function decays as a power law below 



criticality, instead of showing an exponential decay as in 
DP. This is a consequence of having infections that can 
be produced at very large times. We also found that this 
power-law decay is different from the one obtained for 
the two-point correlation function above criticality. 

Subsequently, including ffuctuation effects, we have 
renormalized the theory at one loop. We have calcu- 
lated the renormalization-group equations and we have 
determined the critical exponents at one order in an e ex- 
pansion. The critical exponents vary continuously with 
the Levy parameter and obey at one loop an extra rela- 
tionship with respect to DP, Eq. (|115|l , which is a direct 
consequence of the absence of field renormalization. 

We argue now that the new relationship, Eq. (|115|l . is 
valid to all orders in perturbation theory, which will be 
true if the absence of field renormalization occurs at any 
loop. The absence of field renormalization just by power 
counting is difficult to see in the renormalization scheme 
applied here. This is because terms proportional to E'^ 
in principle can be generated at any loop. Therefore in 
our case it will be necessary to check that the coefficients 
of these terms are not divergent at the upper critical di- 
mension. The absence of singular field renormalization 
to all orders is more clear in other renormalization-group 
schemes — for example, the Wilson method or normaliza- 
tion at nonzero momentum, since in those cases the rele- 
vant Feynman amplitudes are always analytic in the ex- 
ternal energies and momenta. Hence no terms like E'^ as 
E ^ can be generated in loop diagrams, even though 
they are present in the bare propagator. The absence 
of renormalization of such singular terms in the prop- 
agator to all orders has long been known for the case 
of long-range ferromagnets 30], and has also been rec- 
ognized for other variants of DP with long-range spatial 
interactions and with both long-range spatial and 

temporal interactions [sij . 

Finally, we notice that when we set k equal to 1 , we do 
not recover the DP hyperscaling relation from Eq. (|115|l . 
This is due to the fact that the validity of Eq. H115() relies 
on the absence of field renormalization. This does not 
happen in DP, where field renormalization is necessary 
to absorb a divergence that appears at the upper critical 
dimension dc = 4. 

Since this work was completed, a related paper [sj has 
appeared. This differs from the present one in that long- 
range effects in both time and space are considered. The 
renormalization of this theory is simpler than the case 
considered here, because neither the coefficients of E'^ 
nor of k'^ in the bare propagator are renormalized, and 
hence there are two additional scaling relations rather 
than the single one found here. One of the two scaling 
relations obtained in is equivalent to the scaling re- 
lation, Eq. H115f) . found in this work. 
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r(l-K), t ^ oo {p=p^). 

(A5) 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF 
APPENDIX A: ASYMPTOTIC BEHAVIOR OF THE INTEGRALS IN THE RENORMALIZATION 
THE CORRELATION FUNCTION G(^^'(a;,t) CALCULATIONS 



In this appendix we compute how G*^^'^^ (x, t), Eq. (|30|l . 
behaves asymptotically in the limit of t — > oo. The con- 
tribution to the integral is due only to the presence of a 
branch point in i? = 0, since there are no poles in the 
first Riemann sheet. If we call y — \E\, then 



7 — zoo 



1 



dy. 

(Al) 



As t ^ 00, we see from Eq. (|Aip that the leading be- 
havior of the integral comes from the integration domain 
for small y. We consider then a series expansion of the 
integrand as follows: 



1. Integrals involved in the field renormalization 



We proceed to analyze the asymptotic behavior of in- 
tegral / in Eq. (|48|l at the values of u where divergences 
may occur — that is, u — 1/2 and m = ^ ± ioo. Where 
such divergences exist, we will determine the finite coef- 
ficients of the integral. At u ^ 1/2, the numerator and 
denominator of the integrand goes to zero. Expanding 
them around the point u ~ 1/2, we find that they both 
go to zero with the same order in (u — 1/2). Therefore 
the integrand remains finite in the limit u — > 1/2, and 
we infer then that the integral itself does not present a 
divergence in the integration domain around u = 1/2. 

The next step is to consider the asymptotic behavior 
as u — > ^ ± ioo. To analyze this limit it is convenient to 
make a change of variables iy = u — i, and the integration 
/ in Eq. (|49|) can be rewritten as follows: 



1 



1 



2i 



■sin(7rK)y'^ + 0(y2«). (A2) 



(Z?ofc2 + ro)2 
Inserting this result into Eq. (|Aip . we obtain 



'^+'°° dE e 



Et 



7 — iOO 



2Tri E'^ + Dok^ + ro 



t^oo sin(7rK) 1 
7r(i:>ofc2 + ro)2ti+^ 



r(H-K). (A3) 



Therefore, the asymptotic behavior of G*-^'^^ (a;, t) as t 
00 is given by 

r(l + k) sin(7rK) 1 f d'^k e*'^^ 



TT ti+« J (27r)'' (L'ofc2 +ro)2 

-JT^' ^-^00 b<Pc). (A4) 



1=2 I ^y<d/2~2) ( 2 

2tt 



Ay' 



(d/2-2)K/2 



gie((i/2-l)K _ g-ie(£i/2-l)K 



, (Bl) 



where = 9{y) — arctan(2|?/|). We now integrate by 
parts in Eq. IjBip and find 



1 /■°° 
I = - dyy 

71" Jo 



where f{y) is 

fiy)- 



K(d/2-2) 



fiy) = - 

TT 



y 



K,{d/2-2) + l 



dy 



y 



K{d/2 - 2) -f 1 

K(d/2-2) + l 



/(y) 



"K(rf/2-2) + 



Y/'(y) 



, (B2) 



^ ^ (d/2-2)K/2 r ie(d/2-l)K _ ^-i9{d/2-l)K 



4y^ 



(B3) 

The function f(y) is a finite, going to a constant as 
y ^ 00, and vanishes as y ^ quickly enough for the 
first term on the right-hand side to be convergent for 
K{d/2 -2) + l < 0— that is, for d < 4 - - . On the other 
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hand, ior d > 4— ^ , a,s y ^ oo , the first term on the right- 
hand side diverges as y«('i/2-2)+i_ -y^g }ja,ve assumed that 
such a divergence was absorbed into a renormalized mass 
rn in the mass renormalization procedure: the change of 
variables performed in Eq. H47|l scales out of the integral 
the energy dependence in the second term of F^^^^-', and 
therefore the divergences corresponding to mass renor- 
malization in r*-^'^' and Z7^ are the same. Thus, it is 



only the second term on the right-hand side of Eq. ljB2|) 
that may diverge at d ^ dc = 6 — ^. There is already 
a pole at d > 4 — in this term, and thus in order to 
identify the next divergence we will extend the value of 
this term to d > 4 — ^, applying analytic continuation. 

First of all, we need to determine the shape of f'{y) 
for large values of y. We start by considering the series 
expansion of 9 = arctan(2y) = f ~ ^ + 0{y~^) valid for 
y > 1/2. Inserting it in the expression of /(y), we find 



/(y) = 



sin[f(f 



(f-1) 



tan 



tan[f (I - 1)k] 



O 



(B4) 



and 



Ie2 — 



r+ioc 



dE' 
{2m) 



X {d^4)d~ E''^''\E - E'T^i-'^\d^2)d 
{E - E')^{8 -6d + d^)]/[E"' -{E- E'Yf 



E=c 
(B7) 

Let us begin to analyze the existence of divergences in 
Iei- We make a change of variables E' = Eu, where E 
is real and positive. Inserting this change into Eqs. ljB6|) 
and l|B7l) , and after evaluating in the normalization point 
-E = Ci it is possible to see that the only divergences in 
the integrands could appear in the limits u — > i ± ioo. 
In order to analyze these limits, we make the change of 
variables iy — u — ^ and Iei and Ie2 look as follows: 



' d ( 1 ^ 
dyy<-^-^^[x + —^ 

{d - 2){cos[6Ik(| - 1)] - cos[6Ik(| - 3)]} 



An 



sm\0K) 



(B8) 



It is straightforward to see from here that f'{y) oc p- + 
0{^). Then, if we call a = n{d/2 — 2), the integration 
in the right-hand side of Eq. (|B2|I behaves as 



y«+V'(y)d2/ 



oo. 



(B5) 



This result suggests that the next pole would be at a = 0. 
Nevertheless, at the upper critical dimension c?c = 6 — 
2/k, a is different from zero and takes negative values for 
any k, with 1/3 < At < 1. This means that the integral 
in Eq. (|B5|) is convergent as y — > oo. We show with this 
calculation that there is no other divergence present in 
Eq. 



2. Integrals involved in the diffusion constant 
renormalization 



In the Sec. IIVBI the expression of the renormalization 
constant contains the integrals 



Iei 



dE' 



—lOQ 



{2m) 

-E"'{E-E')'''^i-^\d-2) 



2i;'«(5-i) + (£;_£;')«(|-i)(d-4) 
/[E"" -{E- ETf 

E=c 
(B6) 



and 



'B2 



«(4-3) 



8^ ,;o 
X I 2(d-4)dsin 



dyy 
en{ ^-2 



1 



Kf-3) 



4y2 

{d~2)dsm 



(16-6d + fi^)sin(6lK^ ) |./ sin^(6lK), (B9) 



respectively, being 6 = 9{y) ~ arctan(2|j/|). Therefore 
as y ^ oo the integrals diverge as ^^^^^ ^j^g 

upper critical dimension these divergences become loga- 
rithmic. We cannot solve the integrals analytically, and 
for this reason we only determine the coefRcients of such 
divergences. We start by expressing the entire integrand 
as a function of using y = ^^-^ and 1 + = (sin 6')^^. 
For instance, Iei now looks 



Iei = - 



1 



47r 2'=(5-3)+i 
1 

4tt 2«(i-3)+i 



d0f{e) 



d0- 



1 



{d - 2){cos[6'K(f - 1)] - cos[0k{^ - 3)]} 
sm^{9n) 



(cos 61) 

'd 

> 2 



k(4-3)+2 



(BIO) 



The limit of interest is y — > oo or, equivalently, —^ tt/2. 
Therefore, to proceed with our analysis, we can write 
= 7r/2 — a, where a is a variable which tends to zero. 
Substituting as a function of a in Iei in Eq. (|B10|) and 
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taking the limit a — s- 0, we have 



previously, considering E a real and positive number: 



dOfiO) 



da- 



1 



cos(f -a)]«(^-3)+2 



"'0 

id - 2){cos[(§ - aUl - 1)] - cos[(f - aMi ~ 3)]} 
sin^[(^ — q;)k] 
{d - 2){cos[f - 1)] - cos[f - 3)]} 



r-A 

X / da- 
'o c 



A(|-3)+2 



finite, (Bll) 



where < A < 7r/2. We can write k(| - 3) + 2 1 - e/2 
with e ^ dc — d, and then we see that as a — s- 0, the 
integrand / in function of a diverges as 



/(a) ~ aa~ 



a ^ 0, 



(B12) 



with a 



(d-2){cos[fK:(j-l)]-cos[fAc(j-3)]} 



a constant 



number. If we now add and substract this divergence 
in the expression of /(a), we can obtain the coefficient 
of the logarithmic divergence as follows: 



k(|-3) + 1 



du ^"(i-2)-(l-u)"(^-^) 
(27ri) (1 — u)'^ — u'^ 

(B16) 

The integrand is not divergent as w ^ ^, but goes to a 
constant as one can see from a Taylor expansion of the in- 
tegrand. We can infer then that the possible divergences 
may have their origin in the limits it — > ^ ± ioo. For the 
purpose of studying the behavior of the integral in those 
limits, we rewrite the integration in Eq. I)B16(I . making 
the change of variables iy = u — 1/2, as follows: 



du u^(|-2)-(l-7i)'^(i-2) 



1 



(d/2-3)fc/2 
"gie(d/2-2)fc g-ie(d/2-2)K" 



X 



(B17) 



daf{a)= i da[f {a) — aa"'^ ^] + a daa"^ ^ 
Jq Jo 

~a— , e^O. (B13) 



Inserting this result into Eq. (|B10|) . we then have 

-1 C""* 



Iei 



^2«(|-3)+2 KC 



{d - 2){cos[f K(f - 1)] - cos[f K(f - 3)]}' 



sin^dK) 

Preceding in the same way for Ie2 we have 
1 ' 



£ 0. 

(B14) 



E2 



— (c? — 2)c?sin 
dn 



2(d-4)dsin 



TT fd 
2H2-' 



(16-6(i + r)sin k 



—k[ 4 

2 V2 



e 0. 

(B15) 



3. Integrals involved in the composite operator 
renormalization 



In the limit y 00, the integrand diverges as yi^(d/2-3) ^ 
and at the upper critical dimension this integral be- 
come logarithmic divergent. We analyze then this limit 
by making use of the relations y = and setting 

6* = 7r/2 — Of, such that a — * 0. We express Eq. ljB17|l as 
a function of a: 



2 1 ^y«('i/2-3)n ^ J_)(d/2-3)K/2 

2n 4y2 







je{d/2-2)K _ -ie{d/2-2)K. 



-1 



7r2'^(f-3)+i 



sin[f - 2)] 



da- 



1 



,l-Ke/2 



finite. (B18) 



The integrand diverges as a"'/^^^ in the limits of a ^ 
and e — > 0. In order to identify the coefficient of this 
divergence, we add and substract the divergence itself 
from Eq. ljB17|l . finding 



^+*°° du u«(#-2)-(l-w)''(i-2) 



(27ri) (1 - u)« - it« 



a h finite, e -> 0, (B19) 

KC 



The integration in E' in Eq. (|74|l can be studied start- 
ing by doing the change of variable E' = Eu, as we did 



where a is a constant equal to a = rr—. — smf ^ k( 2)] ^ 
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